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METHOD AND APPARATUS FOR IMPROVING RESOLUTION 
IN SPECTROMETERS PROCESSING OUTPUT STEPS FROM 
NON-IDEAL SIGNAL SOURCES 

5 CROSS-REFERENCE TO RELATED APPLICATIONS 

This application claims priority from provisional application 60/151,020 
filed August 27, 1999, the disclosure of which is incorporated by reference in its entirety 
for all purposes. 

10 STATEMENT AS TO RIGHTS TQ INVENTIONS MADE UNDER 

FEDERALLY SPONSORED RESEARCH AND DEVELOPMENT 

The U.S. Government has rights in this invention pursuant to Contract No. 
DE-FG03-97ER82510 awarded by the Department of Energy. 

1 5 BACKGROUND OF THE INVENTION 

The present invention relates generally to signal processing systems, and 
more particularly to processing the step-like output signals generated by non-ideal, 
nominally single-pole ('*N-1P") devices responding to possibly time- varying, pulse-like 
input signals of finite duration, wherein the goal is to recover the integrated areas of the 
20 input signals. 

The specific embodiments described relate to processing step-like signals 
generated by detector systems in response to absorbed radiation or particles and, more 
particularly, to digitally processing such step-like signals in high resolution, high rate 
gamma ray (y-ray) spectrometers with resistive feedback preamplifiers connected to large 
25 volume germanium detectors. The application of measuring the step-like output signals 
from y-ray detector preamplifiers to measure the y-rays' energies is just a specific 
example, and is described because this was the area in which the method was first 
developed. 

The techniques that we have developed solve this problem generally, and 
30 therefore should not be construed as being limited to this specific application. Any 

detection system, for example, that produces output current signals that are integrated by 
charge sensitive preamplifiers could be treated by these techniques, whether the detected 
quantities are light pulses, x-rays, nuclear particles, chemical reactions, or otherwise. The 
techniques, in fact, is not limited to "detector systems"/?^'' are, in fact, general 

35 purpose signal processing techniques which may be broadly applied, once understood. 

The outputs from superconducting bolometers, for example, produce step-like signals that 
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are readily treated by the invention. The field of gamma spectroscopy, where 0.1% or 
less makes the difference between a bad and a good detector, however, provides 
particularly stringent tests of our techniques. 

The term "step-like signal" also requires some discussion. The output of 
5 an ideal single-pole ("IP ") device to an ideal impulse (delta) function input, is an 
infinitely fast rise time followed by an exponential decay whose time constant Xd is 
characteristic of the pole. Viewed on a time scale short compared to xd, this output will 
look like a pure step, while, when viewed on a time scale long compared to Xd, it will look 
like a pulse. A real IP device output, however, will have a finite risetime, Xf, whose 

10 duration will be determined both by the nature of the device and, particularly, by the 

duration of its real input signal. Provided that Xf is significantly shorter than Xd, a real IP 
device output signal, viewed on a time scale comparable to Xd, will then show a risetime 
region, whose shape may be difficult to describe mathematically, followed, after a period 
comparable to Xr, by an exponential decay with time constant Xd- The output of a N- IP 

15 device will be similar, with additional distortions. We will refer to such signals, viewed 
on this time scale, as "step-like". 

Gamma-rav (y-rav) detection requirements 

The detection and measurement of y-ray energies is a well-established 
discipline whose primary goal is to accurately determine both the niunber and energies of 
20 y-rays emitted from some target source. The requirements of good energy resolution and 
high count rate capability usually conflict, however, since count rates are enhanced by 
increasing detector volume, which increases output signal distortion and so degrades 
energy resolution. High count rates also degrade energy resolution directly due to 
practical problems in preamplifier design. 

25 Description of the problems 

The field of y-ray detection is highly developed. A fairly comprehensive 
introduction to the state of the art may be found in the volume "Radiation Detection and 
Measurement, 2nd Ed." by Glenn F. Knoll [KNOLL-1989]. Below we note only the 
issues relevant to the present invention. In the first section, we discuss how pole/zero 
30 cancellation errors introduce a second pole, spoiling the preamplifier's single pole 

response. In the second section, we examine how the finite input signal duration, in this 
case due to charge collection, distorts the preamplifier's output fi'om the ideal. 

Pole/zero cancellation errors 

Fig. 1 A shows a typical solid state y-ray spectrometer comprising a 
35 semiconductor detector diode 7 biased by a voltage supply 8 and connected to a 

preamplifier 10 comprising an amplifier 13 with a feedback capacitor C 15 and resistor R 
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17. As drawn, preamplifier 10 is a single pole circuit whose response to an impulse (delta 
function) input is A exp(-t/x2), where T2 = RC and A is the area under the impulse. 
Because T2 is typically of order 1 ms, which is too long for the following circuits, a 
pole/zero (P/Z) network 20 cancels the pole at 1/t2 and replaces it with a pole at 1/ti, 
5 where xi typically is 50 ^is. Gain stage 22 then amplifies and buffers the preamplifier's 
output signal for shaping amplifier 23 which feeds multichannel analyzer (MCA) 24. 

If the time duration of the current pulse arising from the charge deposited 
in detector 7 by a y-ray absorption is very short compared to xi, the output of stage 22 
will be an exponentially decaying step whose amplitude is the pulse integral and 

10 proportional to the deposited charge, y-ray spectrometers are therefore designed to 
measure these step amplitudes to measure the charge deposited by the absorbed y-ray. 
Other forms of radiation, including neutrons, alpha and beta particles, and x-rays behave 
similarly and their energies are measured the same way. 

Commonly, however, both the input's finite duration and the pole-zero 

15 circuit's imperfections distort the preamplifier's response, destroying the proportionality 
between the output step's amplitude and the deposited charge and so degrading the 
system's energy resolution. Imperfections in P/Z network 20 arise fi-om difficulties in 
precisely canceHng the X2 component, leaving a small residual firaction, of order 1-2%, in 
the output signal. Fig. IB shows a 5% residual X2 component for ease of viewing: an * 

20 exponential decay signal 25 with time constant X2, input to the P/Z network 20, produces 
either output signal 27 or 29, depending upon whether the residual X2 term is positive or 
negative. 

These X2 residuals are particularly bothersome at high coimting rates, 
where each signal step rides upon a X2 background fi-om all preceding steps. As these 
25 arrive randomly, the resulting baseline bias also fluctuates randomly in time, which the 
spectroscopy amplifier's baseline restoration circuit cannot track well. These terms, 
which may only be a few tenths of 1%, become a significant resolution degradation at 1 
MeV where 0.05% energy resolution is desired. 

Signal risetime fluctuations and ballistic deficit 

30 Fig. 2 shows a preamplifier 10 fi"ont end with a cross sectional view of the 

detector 7 of Fig. 1, for the common coaxial geometry. The dashed lines show electric 
field line within the detector body 30, which vary considerably with local geometry. Two 
factors cause charge collection time variations within the detector and thus risetime 
variations in the preamplifier's 10 signal output: 1) the difference between carrier 

35 velocities; and 2) the existence of different path lengths within the detector. 

RAUDORPH-1982 describes these issues. These risetime variations produce ballistic 
deficit by two paths, one direct, one indirect. The direct effect is well understood, per 
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GOULDING-1988: the output filter's response varies with the time dependent shape of 
the charge arrival, being the convolution of the two. A trapezoidal filter greatly reduces 
this effect in the absence of exponential decay. 

The indirect effect source of ballistic deficit is due to fluctuations in charge 
5 loss through the feedback resistor with differing risetime, as seen in Fig. 3A with two 
risetimes, 40 and 42, where Fig. 3B enlarges their peak regions. The slower risetime 
signal loses less charge and thus is larger once charge collection is complete. Even filters 
which ignore the charge collection region are still sensitive to this lost charge effect, and 
relatively small errors of this size can substantially degrade resolution. For a trapezoidal 
10 filter, the collection time difference shown Fig. 3B produces a 0.2% amplitude difference 
(2,000 eV at 1 MeV) which will degrade ideal 1 .7 keV resolution to 2.6 keV, Ballistic 
deficit errors must therefore be reduced to less than 0.05% to obtain ideal spectrometer 
resolutions at 1 MeV, 

Charge Trapping Losses 
15 Charge trapping also produces errors in y-ray energy measurements since 

trapped charges are lost to the measurement. The present invention does not seek to 
address this problem. 

Generalizations 

It is important to note that the pole-zero cancellation errors described 

20 above do not arise firom the preamplifier's use in y-ray spectroscopy, but are a generic 
problem in low noise, charge sensitive preamplifiers. Similarly, while the described 
risetime fluctuations described above were attributed to the geometry and construction of 
large volume Ge detectors, it is clear that such problems fundamentally arise from 
interactions between the finite charge collection time and the electrical characteristics of 

25 the preamplifier and not from the physics of the collection processes. Geometric heat 
flow variations in the photon absorbing mass produce similar effects in the 
superconducting bolometers mentioned earlier. Risetime issues are therefore a potentially 
general problem as well, and may need to be corrected for in other, non-y-ray; detectors 
whenever the highest measurement accuracies are required. The methods we describe 

30 offer just that capability. Further, the terminology "single-pole" or "multiple-pole" 
comes firom the LaPlace Transform treatment of differential equations describing time 
variant phenomena. Any device which shows "single pole" behavior, for example, will 
display exponential time decay in response to an impulse input and, therefore, may be, for 
example, mechanical, thermal, chemical, or magnetic in nature in addition to the 

35 electronic case presented here. Our method can be directly applied to these devices as 
well, as will be apparent firom the teachings herein. 
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Existing Correction Schemes 

The prior art deals with pole/zero errors in two ways: first, by canceling t2 
as accurately as possible; and, second, with baseline restoration schemes which try to 
track the shaping amplifier's "no signal" output as closely as possible, an approach which 
S degrades as rates becomes high. We have not foimd any approaches which measure 
and/or correct for the effect directly. 

Over the years, various heuristic schemes have been developed which 
attempt to compensate for ballistic deficit. RADEKA-1982 introduced trapezoidal 
filtering and developed a time-variant implementation, using a gated filter following a 

10 semi-Gaussian shaper, that provided significant resolution improvements. WHITE- 1988 
proposed a different gated integration approach using a series switch to excise the charge 
collection region out of the preamplifier signal entirely. The final circuit was complex 
and had enhanced deadtime problems. GOlJLDING-1988, RAUDORF-1982, and 
SIMPSON-1990 disclose schemes that depend on directly measuring the signal's risetime 

15 tf and correcting the energy filter output by a term like t^". These approaches are complex 
to implement and require precise expert adjustments to operate. The underlying 
assumptions are not particularly valid and improvements in energy resolution have been 
modest in practice. 

HINSHAW-1991 and KUMAZAWA- 1998 describe attempts to correct for ' 
20 ballistic deficit by capturing peak amplitudes fi*om two filters which respond to the 

ballistic deficit differently, one an energy measuring filter and one a differentiating (or 
bipolar shaping) filter. T)T>ically a significant fi-action of their difference in peak heights 
is added to the energy filter's peak to correct it. 

Related Art 

25 There is some related art wherein the details of the shapes of the 

preamplifier output signals are sampled digitally and used either to distinguish between 
different types of particles absorbed in the detector (e.g., MILLER-1994) or to distinguish 
between single and multiple interaction events in large germanium detectors. See, for 
example, TAKAHASHH994 and AALSETH-1998. 

30 WARBURTON-1997, WARBURTON-1998, and WARBURTON-1999 

describe methods for implementing digital filtering and x-ray spectroscopy. While they 
do not address the issues under consideration, some of their filtering techniques will be 
employed in the present invention and are referenced in the specification. 
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SUMMARY OF THE INVENTION 

The present invention provides techniques for measuring a step-like output 
signal from a nominally single-pole (N-IP) device in response to a pulse-like input signal 
to detemiine the integrated area of said input signal. The invention addresses the 
5 possibility that the device deviates from ideality due to the presence of additional poles, 
zeros and/or a DC offset and that the input pulses may have finite time durations and 
variable amplitudes. In a specific example, the measurement determines collected 
detector charges from step-like preamplifier output signals in the presence of both 
risetime fluctuations and imprecise pole/zero cancellations. 

10 In brief, the present invention contemplates processing the N-IP device's 

(e.g., preamplifier's) output signal using a set of one or more shaping filters. Where 
plural filters are used, they typically have different time constants. A set of samples of 
the outputs of this filter set is captured in such a manner that the multiple sample values 
bear prescribed time relationships to one other. We refer to this set of sample values as a 

15 correlated multiple output sample (or "cMOS" for short). The relationships among the 
individual sample values may be determined by the times of their capture, by delay 
elements inserted in the signal paths, or some combination of the two. Further, the 
different individual sample values can be obtained from associated different filters, or 
plural sample values can be obtained from the same filter, but captured at different times 

20 or with different delays. The term cMOS is intended to cover these multiple possibilities 
and is discussed further in § 5.1 below. 

The input pulse's integrated area is determined by capturing a cMOS in 
response to detecting a step-like signal (sometimes called an event). A weighted simi of 
the individual sample values in the cMOS (sometimes referred to as cMOS values) is then 

25 formed as a measure of the input pulse's integrated area (e.g., total charge). The 

weighting factors can be computed directly from information about the N-IP device's 
decay constants, the filter set, and the prescribed time intervals in the cMOS. 

While embodiments of the invention use triangular and trapezoidal filters, 
the invention does not require specific filter shapes. Further, the invention does not 

30 require that the filter's capture times be precisely located relative to the step-like signal's 
leading edge. Rather, the method derives its accuracy from repeatably reproducing the 
set of prescribed time relationships between the sample values in the cMOS. The 
underlying capture method is therefore time-based rather than amplitude-based, which 
differentiates the invention from prior art methods which capture filter samples based on 

35 their maximum amplitudes. 

While the preferred implementation uses digital signal processing, 
implementations using solely analog processing or hybrid approaches are also feasible. 
Embodiments of the invention require neither a measurement of the step-like signal's 

6 
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risetime nor more precise information about its arrival time than is found by the pileup 
inspection circuits of commonly available commercial shaping amplifiers. 

When using a single cMOS, however, the measurement's noise is 
increased by the number N of additional filter samples required to compensate for the 
N-IP device's non-ideal terms. For those cases where this increased noise is significant, 
we also show how N may be reduced by creating a parameterized model of the non-ideal 
terms; making baseline measurements to determine the parameters; and then using them 
to correct area (e.g., charge) measurements made using simpler sets of filters which have 
less noise. In some implementations, these baseline measurements are cMOS values 
captured at times when the filters are not processing the step-like signals. In the detector- 
preamplifier case described in detail, the resultant implementation achieves the low 
electronic noise levels of conventional trapezoidal filters while also eliminating resolution 
loss due to both ballistic deficit and high count rates. 

A further understanding of the nature and advantages of the present 
invention may be realized by reference to the remaining portions of the specification and 
the drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 A shows a schematic of a semiconductor, y-ray detector fi-ont end, 
including a resistive feedback integrator, pole/zero network and output gain stage; 

Fig. IB shows voltages generated by the preamplifier portion of the circuit 
in Fig. 1 A when the pole/zero circuit is not ideally adjusted; 

Fig. 2 expands the detector structure of Fig. lA to demonstrate the 
geometric and other factors which contribute to signal risetime variations; 

Fig. 3A shows the variation in risetime between pure electron and pure 
hole collection for a simple collection path in the detector structure of Fig. 2; 

Fig. 3B shows an expanded view of the signal peaks in Fig. 3A; 

Fig. 4 shows a step-like feature in a general preamplifier output, broken 
into its component signals, and establishes a nomenclature system for describing these 
signals mathematically; 

Fig. 5A and 5B show how to compute the appropriate intervals between 
capture times of the multiple filters; 

Fig. 6A and 6B show a pair of purely digital circuits which implement the 
described method; 

Fig. 7 shows a purely analog circuit which implements the described 

method; and 
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Fig. 8 shows a hybrid analog/digital circuit which implements the 
described method. 

DESCRIPTION OF SPECIFIC EMBODIMENTS 

1. Introduction 

5 The invention has three basic steps: applying a set of filters to the 

preamplifier output; capturing a correlated multiple output sample ("cMOS") fi-om the 
filter set in response to a detected event; and forming a weighted sum of these sample 
values to accurately recover the energy of the detected event. The cMOS is a set of 
samples bearing a prescribed set of time relationships between one another. The 

10 appropriate time relationships are achieved by setting the times of capture, inserting 

delays into the signal paths, or a combination thereof, and are independent of the specific 
risetime shape of any particular event. Our overall goal, therefore, is to describe how to 
select an appropriate set of filters; explain what it means to capture a cMOS and how to 
achieve it; and show how to compute the weighting functions. 

15 The theory underlying the invention method is based on modeling the 

preamplifier's (or N-IP device's) step-like output response (the "output") to the input 
charge (or impulse: the "input") generated by the event as the sum of a small number of 
analytical functions which are represented by the same number of amplitudes. In a 
specific embodiment, five amplitudes are used to describe: the xi and T2 exponential 

20 components prior to the step-like signal, the ideal amplitude of the step itself, an error 
term due to the balUstic deficit, and a DC level. We first reduce the number of 
amplitudes to four by making differential filter measurements which are DC level 
insensitive. To find these four amplitudes, we make four independent measurements 
(fi"om three filters) in the step's vicinity, express each measurement as a linear 

25 combination of the four amplitudes, and solve the resultant set of four linear equations to 
recover the ideal step amplitude, which is proportional to the input charge (or impulse 
area). Obtaining a set of equations which is well conditioned and insensitive to 
preamplifier noise, ballistic deficit, and T2 amplitude is then primarily a matter of 
selecting an appropriate set of filters and time correlated capture points, which then define 

30 the CMOS. 

Step amplitudes obtained using the four-filter measurement method have 
almost a factor of 3 higher electronic noise than simple trapezoidal filters. In cases where 
this is important, we show how simpler filters can be corrected for the 12 and DC 
components by developing a parameterized model of these terms and making one or more 
35 baseline measurements to determine their values. When multiple measurements are 
made, the parameters must be compensated for their evolution in time, creating a time- 
compensated model. 

8 
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The description below is organized as follows In § 2, using a static, 
captured step-like signal, we locate several filters in the vicinity of its leading edge and 
mathematically describe their outputs in terms of the base amplitudes. In § 3, we use 
these results to set up and solve the set of linear equations to obtain the input charge. 
5 This solution has the fomi of a weighted sum of the filter values, showing the genesis of 
the weighting coefficients. In § 4, we evaluate the method for two filter sets which 
compensate simultaneously for inexact pole/zero cancellation and ballistic deficit. In § 5, 
we introduce the concept of the cMOS to show how to implement the invention in real 
time processing situations. In § 6, we will describe various circuits which implement this 
10 method: a purely digital pair, a purely analog circuit and a hybrid analog-digital circuit. 
In § 7, we analyze our method's electronic noise. Finally, in § 8, we show how a simpler, 
three filter measurement can be used to achieve similar accuracy with lower noise, 
provided time recorded baseline measurements are made to correct for the T2 and DC 
preamplifier components 

15 2. Mathematical description of generalized filter outputs 
2. 1 . Decav components of the step-like signal 

Fig, 4 shows a typical preamplifier step-like output signal as it might look 
if it were digitized by an analog to digital converter (ADC) and saved in a computer 
memory. The shown time scale is arbitrary, with zero set a little before the onset of 

20 charge collection for convenience. To simplify the image, only a few members of the set 
of discrete values have been identified explicitly: subscripted capital Q's denote values at 
specific instants since the signal really represents the charge integrated on feedback 
capacitor 15 by amplifier 13. The T] component has three regions: an exponential decay 
region 50 prior to the step (i.e., times less than 0), a charge collection region 52, and an 

25 exponential decay region 53 following the step (i.e., times greater than 1). The curve 55, 
showing how curve 50 would have continued to decay had the step not arrived, represents 
the continued time decay of all previous steps, that is, the preamplifier's first pole residual 
response to those previous steps. The signal's T2 component has the equivalent four 
regions 57, 58, 60, and 62. Finally, there is the preamplifier's DC offset value Bq 65, 

30 which is constant. The xi and T2 charges injected by the y-ray pulse decay independently, 
since the preamplifier is a linear device. The total collected charge Qgj generating this 
event thus divides into Qigj = (l-P)QgT for xi and QzgT = PQgT for X2, P being the 
division ratio. 

It is important to note that, in the following, we will be describing a signal 
35 processing technique which relies on the mathematical form of the signals and not on 

their physical source. While we develop the method in the context of preamplifier signals 
fi-om y-ray detector preamplifiers, the method itself can be usefully applied to any signals 
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which have the appropriate form. Similar outputs would generally be obtained from 
nominally single pole devices (N-IP devices) in response to impulse inputs of extended 
duration and variable amplitude. In the present case these signals result from a single 
pole electronic circuit responding to an input pulse from a detector, in other cases the 
5 detector itself may produce signals of the appropriate form. The temperature response of 
a superconducting bolometer, for example, shows this behavior [STAHLE-1999] and 
could be treated by the method. 

2.2. Running-sum filter output 

If r values are sampled at uniform intervals At from an exponential decay 
10 with time constant xj, then the value Qjj is given by: 

Qj,i = Qj,oexp(-i At/Tj) = Qj,obji,for(0 i r-1) (1) 

where bj = exp(-At/Tj) is the decay per At. Therefore, a digital filter, defined by 
weighting constants {wj^i^, 1 i r-1 }, acting on the values in Eqn. 1 gives: 

a jk,r = Qi, }Vijc = Q io b >i,k = Q jo A*jk^ 

15 where the * marks this as a general filter result. This shows that the filter output is linear 
in the amplitude Qj o, multiplied A*j^k,r> which is a constant for a particular filter and 
decay constant xi. For a simple miming-sum filter (see WARBURTON-1999), w equals 
unity, so 

l-b^ 



^>=ll'''i=T^;-i(>-'>'i) 



(3) 

20 Eqn. 3 also defines the division factor rj, which is equal to the sum to infinity. 

2.3. General trapezoidal filter output with risetime in its gap 

Next, consider a running-sum oj^^ of m samples starting at time 1 over the 
signal in Fig. 4. Decomposing the signal into 5 components with amplitudes {Qgi3, Qis, 
Qg23, Q23» and Bo}, Eqns. 2 and 3 give: 

25 = QgnA 1 ,ni + Q 13A 1 ,n, + Qg23A2,m + QlZ^Xm + ^ 

where the index j in Aj^m equals 1 or 2 for xj and T2. The equivalent sum 0£-, starting at 
time -4, which only sees amplitudes { Q13, Q23, and Bq} is then: 

= QkAi^ + Q20A2.m ^0 . (5) 

The difference (ct£-h - as.) is a trapezoidal filter of risetime m and gap g. Since Qio = 
30 br"»-g Qi 3 and Q20 = h2'^-^ Q23 : 

= - ^r- = QgnA i,ni + Qi3a(l,m,-m-g) + Qg23A2^ + Q23a(2,m,-m-g) ^ 
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where a(j,r,s) = Aj,r(l-bj^) and the DC components cancel. In Eqn. 6, we stress that Q13 
represents, at time instant 3, the preamplifier's dominant pole (xi) response to all 
preceding step-like signal events; Q23 similarly represents the secondary pole (T2); while 
Qgl3 and Qq23 are their amplitude responses to the present event. Eqn. 6 takes its simple 
5 form because the signal rise occurs entirely within the filter gap region g, so only 

decaying exponential terms are measured outside. Thus as is relatively insensitive to the 
precise location of the step within the gap region because its weighting constants are 
negligible (zero, in fact) there. 

2.4. Triangular filter over the charge collection region 

10 To explicitly probe the step's risetime region, we will use a triangular filter 

(gap = 0) Qg = (ag+ - ag.) of risetime g, whose subtractive leg ag. sits precisely within the 
gap of ay. (This relationship will also be easy to achieve when we implement a digital 
real time system.) We will show that the precise location of the step within the gap is not 
important. Then, as above: 

= QgiaA 1^ + Q j3a(l,g,.g) + Qg23A2,g + Q23«(2,g,-g) - Z q ,gj - Z q2g,i 
1^ ' ' , (7) 

where the qjg^j are the sampled charge values at each interval in the subtracted running- 
sum region, each comprising the new charge increment collected in that interval plus the 
exponentially decaying remnants of the charges collected in the previous intervals. 

Since each charge increment decays to zero independently, while their sum 
20 equals the total collected charge, we can write 

r,Q,gT=.Eqig,i ^g^^ 

00 

'•2Q2gT=.i;q2g.i ^^^^ 

Splitting the infinite summation, however, also gives: 

ri QigT= r, QgT(l - P) = q ig.i + r, Qg,3 

25 so that 

J-i QigT= r , QgT(l - P) = %^ig,i + n Qgl3 



(9) 



(10) 



The correction charge Qc defined in Eqn. 10 is the difference between the total x\ 
collected charge and the amount remaining at the end of gap g. Qc will typically be of 
order 2% of Qgj. As a small correction term, Qc will be relatively imcorrelated to Qgi3, 
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SO the set of linear equations we solve for Qgis and Qc be well conditioned under 
inversion. 

Similarly, for the X2 component: 

2, q2g.i = r^hgT P - Qg23} - r2 Q,2 

5 where Qc2 will be smaller than Qc by approximately the ratio xi/t2. Eqns. 10 and 11 
allow us to express Qg entirely in terms of Qgj, the total charge, and values from outside 
of the charge collection region, where the signals behave strictly as decaying 
exponentials. Equivalent integral expressions may similarly be derived in the continuous 
signal case or for filters with other weighting constants. 

10 2.5. Eliminating 

When the gap region is much smaller than both xi and T2, Qc2 becomes 
proportional to Qc to a high degree of accuracy. This may be seen by defining Qgi as the 
total y-ray charge collected at time interval i in the collection region. Since each Qgi 
decays independently, we can write for the two components: 

15 xl(i-P)Qgii>r=Qgi3 



and (12a) 

,lipQ.br.Q,._ ^^^^ 



while, by definition, 



(13b) 



j2.„(l-P)Qgibr=Qgi3 ^ 

and (13a) 

.||,PQgibr=Qg23 
20 Thus, for Qc: 

Qc = ( 1 - P)QgT - Qg,3 = (1 - P).|; Qgid - bf) ^ 

Expanding biS-> to good accuracy as (l-(g-i)At/Ti) gives (l-big-») equal to (g-i)At/Ti so: 

Q, = (l-P)f .|:'Qgj(g-i) 

^•• = 0 ^ (15) 

Similarly: 

Qc2=P^.SQgi(g-i) 
where p2 is equal to At/T2. Comparing Eqns. 15 and 16: 
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Qc2 is therefore quite small, being about 0. 1% of Qc or 0.002% of Qgj- We can therefore 
express any of our filters in terms of the four quantities Qgj, Qc, Qgis and Qg23 using: 

g-i 

i?o^'«.'^'-'Q<= (18a) 
Qgi3=0-P)QgT-Qc (18b) 



Qg23=PQgT-FQc. (18d) 

3. Solving for collected charge Ogx from the captured waveform 
10 3.1. Filter expressions in terms of the base amplitudes 

§ 2.5 shows how all difference filters in the vicinity of the step's leading 
edge can be expressed as linear equations in the four amplitudes QgT» Qc» Qgl3 Qg23 
(non-difference filters will also require Bq). Therefore we must make four (or five) 
independent filter measurements to solve for QgT- While a wide variety of filters could 

15 be used, certain sets will be better conditioned than others and will also require less 
measurement precision to accommodate step risetime variations. Research and 
experimentation have shown the following set to function well: 1) the trapezoidal 
"energy" filter of risetime m fi"om Eqn. 6; 2) the triangular "risetime probe" filter of Eqn. 
7; and, 3) a pair of triangular "xi and T2 background" filters of risetime m/2, situated at 

20 times -4 and +1 in Fig. 4, respectively. It is important to notice that, in this selection, 

only the risetime probe filter processes the step's risetime, which lies in the energy filter's 
gap and is entirely outside the range of the two other triangular filters. This placement, 
where only the risetime probe filter has non-negligible weighting constants in the risetime 
region, assures that their outputs, and hence Qgx estimate, will be insensitive to the 

25 precise location of the signal step within the gap, which, in real time implementations, 
will remove any need for great timing precision. That the filters share some common 
starting and stopping points is primarily a matter of mathematical and implementational 
convenience, but is not fiindamental to the technique. 

Substituting Eqns. 18 into Eqns. 4 through 7 gives four linear equations 

30 for a^H-, a^-, as, and ag in terms of QgT, Qc. Ql3 and Q23: 



= QgT[(l-P)A,,ni+PA2,. 



Qc 



Ql3A,J 
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+ Q23[A2^]+mBo 

^Z- = QgT{0] - Qc[ql + Qi3[b';""^A , ^] + Q23[br~^A2.„,] + mBo 
= QgT[(l-P)A,.m+PA2.m| - Qe[A,.„+FA2^] + Q,3[a(l,m,-in-gil 



(19a) 
(19b) 



— =QgT 

r, 



+ Q23[a(2,m,-m-g)] 
(1-P)A,^+PA2, 



(19c) 



Q23 



(2-b1)+-F(2-b!) 



a(2,g,-g) 



fl3 



a(l,g,-g) 



(19d) 



We divided CTg by ri so that the coefficient of Qc will approximately equal unity for 
inversion stability. 

Using the same methods for the outputs of three filters of length m/2 (two 
triangular filters 043 and 021 whose running-sum first points are: 043 = [time +3: time 
+1] and a2i = [time -2: time -4] and a trapezoidal filter 032 = [time +2: time -2]) gives: 



a43=-Q 



gT 



-Q.3[a('.?.f)]-Q23[4f.?)] 



«^21=-Ql3b 



•l((i-P)A,^+PA2^^ -Q,[a,^+FA2^] 



(20a) 
(20b) 



(20c) 



These filters also exclude the step's rising edge for location insensitivity. 
While this is not a requirement of our invention, it significantly simplifies both the 
coefficient derivations and trigger system timing requirements. 

3.2. The four-filter solution for 

Collecting results for the four filters listed at the start of § 3.1, we have: 
= Qgijc i-P)A , ^+PA2^] - q|a, ^+FA2,„ 



14 



wo 01/17103 



PCT/USOO/23333 



+ Qi3[a(l.m,-m-g)] +Q23[a(2,in,-m-g)] 
(l-p)A,.g+PA2,g 



(21a) 



+ Q 



13 



a(l.g,-g) 



-Q. 



+ Q23 



(2-M)+Jl^(2-b|) 



a(2,g,-g) 



*43 



-Qi3[a('.?.f)]-Q23[4.f.?)] 
^21 =- Q,3[br"" "(^'T'?)] - Q23[b2'^ "(2.^.^1 



(21b) 



(21c) 



(21d) 



as four linear equations in four unknowns. The unknown parameters are of three types: 
1) (Pi) the unknown charge Qgx (the area of the input impulse in the general case); 2) 
(P2) charges describing the step's finite arrival time (Qg) and the residual charges in the 
10 Ti and X2 decay modes (Q13 and Q23) fix>m previous events; and 3) (P3) a set of fixed 
parameters {xi, X2, P, m and g} which describe the preamplifier's (or N-IP device's) 
transfer fimction and the applied filters. Since the coefficients in [ ]'s are functions only 
of the P3 parameters, they are constants which need be computed only once for a given 
preamplifier and filter set. 

15 Eqns. 21 can be solved by several methods, including casting them into 

matrix notation and then inverting the coefficient matrix. Thus: 



*43 



^21 
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10 



15 



20 



(l-P)A,g+pA2^ 



-(2-b^)-^(2-bf) 



a(l4n,-n>g) 

"(1.8.-8) 
r. 



^1 



a(2^,-in-g) 

«(2&-8) 
r, 



'1 



-(1-P)a(i .^)-pa(2,!5.f) a(i .f .f)+Fa(2.^,f) -a(i.f .f) -a(2.f -f) 

0 0 -b7'^^a(i.f.f)-b-'^8a^2.f.5) 



QgT 

Qc 
Ql3 
Q23 



becomes 
so that 

and 



o =J»Q 



QgT--JTl <^r'*" JT2— "^JT 



'13 ^43 ^•'u^ai , 



(22) 
(23a) 
(23b) 

(23c) 



As a solution for Qgj, Eqn. 23c has the advertized form: a weighted sum 
of the four filter outputs. The weighting coefficients Jij-^ are computed only once for a 
particular set of filters and P3 parameters {xi, 12, P, m and g} and then multiplied by the 
four filter outputs as, ag/ri, 043 and a2i to obtain the total pulse charge Qgj (and hence 
the energy of the absorbed y-ray) for each captured event. We note that if more filters are 
measured than there are parameters Q, the matrix can be inverted using least squares 
methods to obtain Eqn. 23b. 

4. Evaluation for two filter sets 

4. 1 . A four>filter set with tvpical parameters 

We now evaluate a typical case of the filter set leading to Eqns, 22 and 23: 
we set m equal to 4 fxs, g to 1 ^s and assume typical values of xi = 50 ^is and X2 = 1 ms 
with a division ratio p = 0,04. In this case, for 40 MHz sampling, m = 1 60, g = 40 and we 
can compute 



J = 



so 
and 



1.54041E+2 -1.54138E+2 -8.08795E+1 -8.0041 lE+1 
1.98087E-2 -1.06150E+0 -2.00007E-3 -1.99953E-3 
-2.95906E+0 3.07604E+0 -1. 53785 E+1 -1.59682E+1 
O.OOOOOE+0 O.OOOOOE+0 -1.69959E+1 -1.60483E+1 



6.309486E-3 -9.623578E-1 -1.593231E-2 -1.549595E-2 
1.190638E-^ -9.600143E-1 -2.284270E-4 -2.469335E-4 
1.194301E^2 6.874687E-2 6.221807E-1 -6.786521E-1 
-1.264822E-2 -7.280624E-2 -6.589193E-1 6.564133E-1 



(24a) 



(24b) 



QgT = 6. 309486E-3 - 4.8 1 1 789E-4 a^-l .59323 1 E-2 a^^ - 1 .549595E-2 , ^24^) 
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cTg is typically of order 0.25a5:, so the CTg term is typically less than 2% of QgT- In spite 
of its small size, however, Qg carries the burden of the ballistic deficit correction, since it 
is the only filter that probes the step's risetime directly. 

The coefficients in Eqn. 24c also depend upon the value p. When the 
5 pole-zero is perfectly balanced, (p = 0) Eqn. 24c becomes: 

Qgj = 6.309486E-3 a^, - 4.81 1789E-^ - 1.593231E-2 cj^^ - 1.549595E-2 o^^ ^25) 

where ag's coefficient has changed by 4%, a^'s by 0.04%, and the a43's and 021 's by 
about 0.008%. The need for p accuracy is therefore a function of how accurately one 
wishes to correct for risetime fluctuations and pole-zero inaccuracies. The J"^ weighting 
10 coefficients are found to vary linearly with p, which allows p to be determined 

experimentally if it is not precisely known. To do this, first select a set of p values 

bracketing its expected value and compute Eqns. 24 to produce a set of weights 
(the coefficients in Eqn. 24C) for each p value. Next collect data for each weight set 

{Jin measure the standard deviation aQgT(P) in QgT (i.e., the energy resolution) 

15 which is obtained. Then, by plotting or interpolating the crQg-KP) results versus P, select 
a P value which minimizes CTQgT(P)- If the preamplifier is stable, p should only need to 
be determined once. 

We tested the accuracy of Eqn. 24c's in the presence of both risetime 
variations and a sizable T2 decay term (equal to the xi step amplitude) using a spreadsheet 
20 to compute the filter responses for modeled wavefomis. Table 1 presents these results. 
As may be seen, the Eqn. 24c is insensitive to both risetime fluctuations and 12 
components at the few ppm level and is three orders of magnitude better than a simple 
trapezoidal filter. 

Table 1 : Modeled charge Qgj for 5 filters in the presence of varying risetimes and X2 
25 amplitudes for the case of p equals 0.02. 

Risetime = 450 ns Risetime = 950 ns 



Model 



<T2> = 0 



<X2> = 1 



<T2> = 0 



<T2> = 1 



30 



Triangular 
Trapezoidal 



0.861966 

(100.00%) 

0.867397 
(100.00%) 



0.859787 
(99.75%) 

0.862421 
(99.43%) 



0.831149 
(96.43%) 

0.863901 
(99.59%) 



0.827170 
(95.96%) 

0.858928 
(99.02%) 
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Section 4.1 
(Eqn. 24c) 

Section 4.2 
(Eqn. 33) 



0.999997 
(100.000%) 

1.000119 
(100.000%) 



0.999997 
(100.000%) 

1.000113 
(99.999%) 



1.000001 
(100.000%) 

1.000028 
(99.991%) 



1.000001 
(100.000%) 

1.000022 
(99.990%) 



4.2. A three-filter simplification which still has ^ood performance 

For cases which do not demand the precision of Eqn. 24c, or have small T2 
components, we show a simpler, three-filter set which still has good performance. We 
use trapezoidal/triangular filters for their insensitivity to the *'DC-like" slow variation in 
the T2 component: namely the same "energy" and *'risetime probe" filters as in § 4.1 and a 
second trapezoidal filter 032 (Eqn. 20c), aligned to exclude the risetime region, to 
estimate Q13. Thus: 



















Zi 




0 


1 


Wa,2 




WA, m 


-WA, m 
'•■3 



-a(l,m.-m-g) 
g a(l>g,-g) 



Wa 



(1 m J 
J '2— -9 



QgT 
Ql3 



(26) 

In Eqn, 26 we reversed the sign of as in order to create a filter that is difTerential on both 
sides of the step's rising edge and we weighted 032 by W to determine the value that 
gives the best results. We find the lowest noise is for W = 1, so we can repeat Eqns. 24a 
and 24b to obtain: 

QgT = 9.338925&-3 a^.- 5,000000E-4 Og H-3.106015E-2 (27) 

Table 1 shows that this result is sensitive to the X2 component at the few 
ppm level and to risetime variations at the 100 ppm level. For many practical cases, this 
is sufficiently accurate. Assuming a random distribution of risetimes, for example, this 
would add less than 100 eV in quadrature to the energy resolution at 1 MeV, which would 
be scarcely perceptible. 

5, Correl ated multiple output sampling in real time implementations 
5.1. Defining "correlated multiple output sampling'' 

The filter weighting Eqns. 25 and 27 derivations were developed as if the 
signals had been captured into memory, so the running-siuns could simply be taken over 
the desired sample values. This implementation is feasible at lower data rates, but in the 
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preferred implementation, as taught by WARBURTON-1997, a gate array is used to 
realize the difference of running-sum filters, which run in real time, each outputting a new 
value its own fixed number of clock cycles after a new ADC sample is supplied. A priori, 
then, there may be no single instant when a set of filter values can be captured which will 
5 reproduce the inter-filter time relationships required by Eqns. 25 and 27. Instead, in 
general, the outputs of the different filters will have to be captured at different times and 
these times will have to be carefully correlated if Eqns. 25 and 27 are to produce accurate 
values of QgT- Because multiple filter values must be captured, we call the resultant set 
of captured values a "correlated multiple output sample" or "cMOS'* and refer to these 

10 values either as "cMOS values" or as "individual sample values in the cMOS". 

We also observe that the number of cMOS values required by Eqns. 25 
and 27 and number of realized filters needed to produce them need not be the same. For 
example, in Eqn. 25, the 021 and 042 filter values both come fh>m triangular filters of 
length m/2 with their starting points separated by the interval m+g. Thus two outputs 

1 5 from the same realized filter can be captured to supply the two required values. 

While we have noted that the filter ranges do not have to be precisely 
positioned relative to the step 's rising edge, the accuracy of their positions relative to one 
other in the cMOS capture determines the accuracy that the method can attain. This is 
because the running-sum ranges, which translate directly into filter lengths and time 

20 intervals between starting locations, are encoded into the J-^ weighting coefficients 

through terms like aO',r,s) = Aj,r(l-bjS) (see, e.g., Eqn. 6), where r is the filter length and 
s the interval between the sum start points. For simplicity, we aligned the various filters - 
so their starting points coincided, but other cases are trivially treated by scaling Eqns. 4-7 
by terms like bj^ where t is the offset between starting locations. Therefore, there is an 

25 intimate relationship between the weighting coefficients and the filter ranges which, if 
violated, will introduce inaccuracies. The 's are computed for a particular set of filter 
sum ranges. If another set of ranges is captured, then expected cancellations between 
filter terms will not occur exactly and the differences will appear as errors in the detected 
charge. In digital implementations, the fact that the signal samples are discrete and 

30 enumerated makes is relatively easy to maintain accurate time relationships by measuring 
time in 'ticks" of the digitization clock and then simply counting out the required time 
intervals. 

The requirement of appropriate time correlations between the cMOS 
values clearly differentiates the present method from prior art as represented by the 
35 approaches of HINSHAW-1991 and KUMAZAWA-1988. While two filter outputs were 
captured and weighted in those methods, the capturing steps used peak detecting means. 
The methods are therefore amplitude based rather than time based, as in the present 
method, and the captured values carry different information about the step shapes. 
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Further, the method of capturing based on peak detection assures that there will be no 
fixed time relationship between the captured samples since the two filters' times of 
maximum amplitude depend differently on step shape, which will vary from step to step. 

In light of the above, it is clear that there will be a wide variety of possible 
approaches to capturing filter cMOS values with the prescribed time relationships. The 
outputs of single filters can be captured multiple times. Filter outputs can be captured as 
they appear from the gate array or can be delayed for convenience. We will show several 
examples in § 6. Which method is best will depend upon implementation details and will 
vary from case to case. We therefore intend the phrase "cMOS" (for "correlated multiple 
output sample") to cover the multiplicity of possibilities for capturing, in a real-time 
implementation, filter values whose time relationships are correct for the set of J*^ 
weighting coefficients that will be applied to them to recover the input charge from the 
detected event. When the cMOS capture is made in response to the detection of a step in 
the preamplifier's output signal, with the intention of using it to determine the charge 
deposited in the detector (in the general case the impulse area) and hence the gamma- 
ray's energy, we will often refer to it as an "energy cMOS capture" (or "charge cMOS 
capture" or "area cMOS capture" as appropriate). 

5.2. Timing considerations 

Filters, in real time implementations, require finite times to compute. For 
example, consider a trapezoidal filter of length (risetime) m and gap g over a set of values 
{qi}, where qn is the last sample value ("end sample'*) to enter the filter: 

n n-nv-g 
F=n-ni i=n-2m-g ^ (28) 

If implemented in gate arrays according to the method of WARBURTON-1997, the 
output <yj:(n) will appear 3 filter clock cycles after sample qn is presented. This is its 
propagation delay. For a 4 |lis peaking time (filter clock equals 200 ns), the propagation 
delay is 600 ns. These delays must be included in the design to assure that the cMOS 
capture constraints are satisfied. 

Figs. 5A and 5B show how to compute sampling times for the Eqn. 24c 
filter set. The calculation has two steps: finding the time delays between their end 
samples of the different filters and then correcting for their different propagation delays. 
For the first step we use Fig. 5A, which shows the filter regions relative to the 
preamplifier output step 72. a£ 74, 043 76, 021 77, and Og 79 were all described earlier, 
at 81 is a fast channel timing filter, as described by WARBURTON-1997, which is used 
to detect step-like signals in the preamplifier's output. The timing filter's output triggers 
a digital discriminator when it is located approximately as shown with respect to the 
leading edge of the signal step 72. We can now measure end sample delays relative to at 
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81 as shown, with the filters as 74, CT43 76, and Og 79 appearing at times At^, At43, and 
Ats after at 81 and a2i 77 appearing At2i earlier. From Fig. 5A, if at has a peaking time 
of T, then At^ equals (m+g-2T). 

It is important to note that, while the firing of at may jitter somewhat 
depending upon the step's risetime shape, this only causes the step location to jitter within 
the gap g. The relative timings of all the filters to at, and thus to each other, are fixed by 
construction. Thus, provided Og is wide enough to catch the step for any reasonable 
discriminator triggering, the filter will work as designed. This relative insensitivity to the 
timing accuracy with which the step-like signals are detected is a significant advantage of 
the present invention. 

Fig. 5B shows the second step: correcting for propagation delays. We note 
that these delays may exceed a microsecond in long time constant filters and typically 
increase with filter length. Fig. SB shows propagation delays with respect to the ADC 
sample time ts 83 for a 40 MHz implementation with g = 1 ^is and m = 4 us, so 6tt = 75 
ns, 5tg = 5t2i = 6t43 = 300 ns, and 6ts = 1,200 ns. Combining Figs. 5A and SB, we 
should sample the output of a given filter ax, at 

tx = At,+St,-5t, ^29) 

following the discriminator firing on at's output to assure the cMOS temporal 
relationships required to extract Qgj. If tx is positive, then a counter can be used to 
measure tx and capture the filter value. If tx is negative, however, then its propagation 
delay 6tx will have to be increased by at least (tx| to meet the demands of causality 
(capturing must follow step detection in time). In digital implementations this is a simple 
operation using a FIFO memory. In our present example, no additional delay is required 
since a typical timing filter peaking time is 100 ns, so At2i = -200 ns and 

t2i = At2i + 6t2i + 5t, = -200 + 300 - 75 = 25 ns, (30) 

which is positive. 

6. Circui ts implementing the cMOS method 

6. 1 . Fullv di gital filter implementations 

Two general approaches are possible: 1) asynchronously capturing values 
as they appear fi-om the sub-filters, per the discussion of Sect 5.2, and 2) extending the 
filters' propagation delays so that they can be synchronously captured. In the latter case 
all pileup inspection is carried out before the capture signal is given, while in the former 
case the values are captured as they appear and the pileup inspection is used to decide 
whether anything should be done with them afterwards. We will demonstrate one 
implementation of each to show the principles. 
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6.1.1. Asynchronous capture implementation 

Fig. 6A shows a preferred implementation of Eqn. 24c. The digital 
processor 90 comprises an ADC 92, a real-time digital processing unit (RT-DPU) 93, and 
a digital signal processor (DSP) 95. Per the method taught by WARBURTON-1997, the 
5 RT-DPU processes data at the ADC output rate while the DSP processes data at the event 
rate (impulse signal rate). The RT-DPU produces values a^, 043, a2i, and ag, while the 
DSP multiplies them by the coefficients J-* to obtain Qgj. 

The RT-DPU 93 also has two parts, a fast channel 97 and a number of 
slow channels 98. The fast channel contains the fast timing filter at 100, a digital 

10 discriminator 102, some pileup test and timer logic 103, a clock 104, and an output 

register 105. Each slow channel comprises a triangular or trapezoidal digital filter and an 
output register with either a digital delay line (if Eqo. 29 is negative) or a timer (if Eqn. 
29 is positive). Notice that 043 and a2u both being of length m/2, are being implemented 
using a single triangular filter 04321 107. 

15 In operation, Ot 100 processes ADC 92 data until a step-like signal arrives. 

Discriminator 102 detects this step and signals the pileup inspector 103, which strobes a 
trigger line 106 and initiates pileup inspection. If the step's arrival time is important, the 
output of clock 104 can be captured to output register 105 by the same strobe signal 106. 
Filter 0432 1's 107 output, delayed by t2i using 21 -delay 108 if required by Eqn. 29 

20 causality, is captured immediately in the 21 -register 109 by this trigger strobe, which also 
starts the 43-timer 110, S-timer 111, and g-timer 1 12 to time t43, t£ and tg. Since tg is 
shortest, the g-timer 112 times out first, capturing filter ag's 115 output in the g-register 
116. Similarly, the 43-timer 110 and S-timer 111 time out at t43 and tv to capture filter 
<^432\ 107 and 117 outputs in the 43-register 118 and Z-register 119. If this step is not 

25 piled up, then all four filter values are ready to be read into the DSP 95 via data bus 121. 
The pileup inspector can be implemented as taught by WARBURTON-1999 relative to 
the length of the oz filter and will not be discussed further. 

As per the cMOS discussion in § 5.1, the invention method does 
specifically require that the number of filters equal the number of filter values to be 

30 captured. Here we have captured 04321 twice. 

6.1 .2. Svnchronous capture implementation 

Fig. 6B shows an implementation of Eqn. 27 wherein the three filter 
values are captured just as the value emerges fi^om the slowest (cte) filter. The primary 
physical difference between this implementation and that of Fig. 6 A is that the timers 
35 108, 109, and 110 have been eliminated and a second delay (g-delay 125) has been added 
to compensate for the differences between t^ and tg and t32 fi"om Eqn. 29. All the 
remaining parts are the same and carry the same reference numbers as is Fig. 6A. Notice 
that the trigger line 106 now connects directly to the output registers, rather than through 
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timers as before. For this topology, the functional operation of the pileup test and timer 
#1 103 is essentially identical to that presented by WARBURTON-1999. 

A disadvantage of this circuit, compared to Fig. 6A, is that relatively long 
delays may be required, which may consume excessive gate array resources. Recent 
5 programmable gates, however, are appearing with blocks of dedicated memory, so this 
restriction is becoming less significant. One significant advantage of this design is that, 
because of the synchronous capture, the pileup inspection logic design is much simpler. 

6.2. Fully analog implementation 

While our preferred implementation is fully digital, the method is not 
10 inherently digital per se. In certain cases, particularly for very fast step risetimes and 
very short filter time constants, or to allow a high degree of integration, it may become 
preferable to use an analog implementation. Fig. 7 shows an analog version of Fig. 6A. 
The topologies are identical up to the point where sample and hold circuits replace digital 
registers to capture the filter values, so the operation of the circuit up to this point will be 
1 5 clear to one skilled in the art, given the teaching presented herein. The derivation of the J 
matrix, of course, will have to proceed using weighting functions appropriate for the 
analog filters used. Then analog multipliers (160, etc.) and an analog adder 185 are used 
to implement the J-^ coefficients and addition of Eqn. 24C. We note that, for an 
application with fixed filters, the multiply-and-add fimction is just an op-amp circuit with 
20 appropriate resistor inputs to its summing node. If the pileup tester 143 deems the value 
good, the gate and shape circuit 145 passes out the adder 185 output in a form suitable for 
multichannel analysis. 

6.3. Hvbrid implementations 

Shaping filter characteristics which are difficult to implement digitally 
25 may be required in some situations where the other benefits of digital processing are still 
desired. In these situations hybrid implementations may be considered. Fig. 8 shows one 
such implementation 190. The topology is nearly identical to that of Fig. 6A. The major 
differences are: 1) that the fast channel 191 is implemented in analog circuitry identical to 
that of Fig. 7 (although the digital circuit of Fig. 6A could as easily be used); and 2) that, 
30 in the slow channel section 192, the digital filters a432i 107, as 117, and ag 115 have 
been replaced by analog filters 04321 193, 194, and ag 195 followed by ADCs 196, 
197, and 198. The time behavior of this circuit is the same as was described for the 
asynchronous digital circuit in Fig. 6A. 
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7. Noise considerations 

The Eqns* 24C and 27 filters have higher noise than simple trapezoidal 
filters since they weight the sampled set of data points with values larger than unity. For 
example, by breaking the 2m-i-g filtered region of Eqn. 27 into its 4 major sub-regions of 
length m/2 (ignoring the gap g for this argument) and calling running-sums over these 
regions a\ through cj4, we can recompose Eqn. 27 as: 

Q^r = - 9.338925E-3 H-2.172123E-2 a3-2.172123E-2 02 + 9.338925E-3 a, . 

If the noise standard deviation per sample is a, then the total noise sig(QgT) will be 

0.5 

= 0.299 a 

(32) 



sig^Qgjj = ^160 ^8.722E- 5 a j + 1 60 ^4.71 8E-4 cj^ = 0.299 



10 which can be compared to a trapezoidal filter of the same length: 

,0.5 



siglQJ = |16U|7^| +160|-;^| I -0.112 a 

(33) 

Thus the Eqn. 27 filter will have approximately 2.7 times as much electronic noise as a 
simple trapezoidal filter. The result for the Eqn. 24c filter is similar. For a 4 jxs filter, 
whose electronic noise can be as low as 150 eV, this effect is significant at energies like 
15 100 keV (Fano noise 415 eV) but not at 1 MeV (Fano noise 1,300 eV). And, as peaking 
times shorten to obtain higher count rates, the problem will worsen at higher energies. 
We therefore consider how to obtain the benefits of correcting for both ballistic deficit 
and pole-zero errors without paying the additional filter noise price. 

8. Treating non-ideal terms with baseline measurements 

20 Our approach is to extend the baseline averaging concept of 

WARBURTON-1997 to treat time varying terms as well. In the particular case at hand, 
we would like a way to measure Q23 using baseline measurements in spite of the fact that 
it both exhibits exponential decay and is increased each time an event occurs. 

8.1 . Derivation of an example filter 

25 We begin by developing a simple cMOS filter, composed of three running 

averages, of lengths m, g, and m, which we make contiguous for the practical reasons of 
simplifying the math and minimizing deadtime, but not for any fundamental reason. As 
in previous sections, the length m filters will primarily serve to determine Qgj, while the 
length g filter will be used to eliminate ballistic deficit effects by probing the signal step's 

30 risetime region. Starting with Eqns. 19 A, 19B, and 19D, these filters are: 
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10 



20 



25 



and 
where 



Ani = (l -P)Ai,m+PA2,m. A^ = Ai,m + FA2,m, r = +Fr2 



(34(1) 



In Eqns. 34 we now have 4 classes Pi of parameters: Classes 1 and 2 (Pi, 
P2), as before, which include Qgj, Qc. and Q13; Class 3 (P3), Q23 and Bq, which describe 
the preampHfier's "non-ideal" behavior, Q23 being the amphtude of its second pole's 
response to any previous steps and Bq being the DC offset; and Class 4 (P4): all the 
remaining fixed parameters (e.g., xi, T2, P, m, g) which describe both the filters' and the 
preamplifier's transfer fimctions. We will treat Qgx, Qc> and Q13 as unknowns to be 
recovered from a particular cMOS filter measurement, with the intention of independently 
determining Q23 and Bq fi'om baseline measurements. Eqns. 34 are simple enough to be 
solved by direct elimination of unknowns to get: 



QgT=Aj 



A 



-aJ[] 



LQ23+MB0 



15 where 



l,ni 



(35a) 



(35b) 



L = A 



2^ 



m+g 



1 - 



A. /S, 

-A. '•^ 

"^2 



^24n ^ 



and 



M = m 



Aw 



m 



A — ''1 




m+g 



(35c) 



(35d) 



The constants K, L, and M depend only on Class 4 parameters, which are assumed to be 
time independent for a particular detector-preamplifier pair and set of filter lengths. Eqn. 
35 A is therefore a transform function (the first bracketed term) between the Eqn. 34 
cMOS filter values a£+, a^., and Og- and the desired detector charge Qgj, minus an error 
term f^Qs^ + MBJ A„, * ^ transform's response to the non-ideal (NI) 

preamplifier terms Q23 and Bq. Using Eqn. 35a thus requires that Q23 and Bq to be 
determined at a time appropriately correlated to the cMOS capture (e.g., time 1 in Fig. 4) 
and substituted into the error term to produce an accurate value of QgT- The cMOS 
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capture itself can be physically implemented using any of the approaches described 
above. The subscript q on the bracketed weighted sum of cMOS values indicates that this 
was an "area" cMOS capture of the area filter set to deteraiine deposited charge. 

8,2. Single measureme nt error term correction 

A simple, although not the lowest noise, method to estimate the error term 
is to measure it locally for each cMOS energy capture by making the same Eqn. 34 
cMOS filter capture in close time proximity but without a step-Uke signal present in the 
filters, so that Qgx is zero. Then, from Eqn. 35a, with Qgj equal to zero, we have: 



LQj3 + MBo| =A 
Jb 



T-1 
m 



r ^g- 



(36) 



where the subscript b is for "baseline", called so since, as shown, this "baseline" cMOS 
filter capture produces the transform's baseline response to the preamplifier's non-ideal 
terms. We may also call such a measurement a "null charge" or "null energy" ("null 
area" in the general case) capture since it is made when QgT is zero. Substituting Eqn- 
34 c into Eqn. 33A gives: 



QgT = Aj 



A 



q 



-Ka. 



A 



(37) 



where the subscript q labels the cMOS charge measurement to be corrected. 

This correction is moderately simple and has been found to work 
acceptably in some cases. However, if Q23 is large, as may be true at high output step 
rates, then time decay effects in [LQ23 + MBoJb may still be significant. This can be 
approximately corrected by measuring the time interval 8t = te - tb (in clock ticks At) 
between the energy and baseline captures and correcting by the decay term b2^^ to obtain: 



A 



-b' 



5t 



A 



Jb 



(38) 

In this method, if Bq is tmly a constant and the time interval 8t between the energy and 
baseline captures is fixed, then the error term MBob2*^, will also be a constant and merely 
produce an offset in the spectrum. 

Eqn. 38 has two significant problems: 1) increased noise, due to the noise 
in the single baseline measurement, which is the same size as in the measurement itself 
and degrades resolution by sqrt(2); and 2) it increases deadtime by at least a factor of two 
if the two measurements are independent and so separated by at least 2m+g in time. It is 
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therefore useful to be able to make multiple, independent measurements of Bo and Q23 
and average them in order to reduce their variance. We present a method for doing so 
below. 

8.3. The time compensated baseline model 

5 In making filter derivations throughout this work, we have assumed that 

we have an accurate preamplifier model in terms of its several decay constants and DC 
offset. We therefore take the logical step and use precisely the same model to correlate 
measurements made at different times so they may be averaged. Thus, if we measure 
Ql3> Q23. or Bo at one time, then their values At clock cycles later will be given by 

10 bi^*Qi3, b2^*Q23. and Bo, respectively, unless a step arrives within the interval At. In that 
case, at the end of the step arrival time (e.g., time 1 in Fig. 4), we increment Q13 by (1- 
P)QgT> Q23 by PQgT» and Bo by zero, after which their exponential decays will continue 
as before. For the purposes of baseline averaging, then, Q13, Q23, or Bo become "time- 
compensated" model (or "TC -model") parameters, whose averaged values are <Qi3>, 

15 <Q23>» and <Bo>, respectively. 

If we make measurements of Bo and Q23 from time to time, we can now 
update our estimate of Bo simply by using the method of WARBlJRTON-1997 or some 
other averaging scheme. For Q23 our TC-model now gives: 

Q23(i) = b^^Q23(i-l)^ (39) 

20 so, if <Q23^i-l represents an "averaged" estimate of Q23(i-1), last measured at time tj-i, 
we update <Q23>i, using a new baseline measurement of Q23(i) at time tj, by the 
following extension of Warburton and Hubbard's Eqn, 7, where Ati = tj - ti-i: 

Q23(i) 

(40) 

Equivalently, we update <Q23> for the arrival of an impulse signal of charge QgiCO at 
25 time ti using: 

(Q23)i = b^''(Q23>(i-,)+PQgT(i) ■ (41) 

If required, Q13 can be similarly extended from its last time of measurement by replacing 
b2 withbi in Eqn. 39. 

<Q23> 2md <Bo> are used with an charge cMOS capture to determine Qgx 
30 using Eqn- 35a as follows. First, <Bo> simply replaces Bo in the equation. <Q23>, 

however, has to be extended from its time of last update ti-i to the present capture time ti 
using Eqn. 39, according to our TC-model, where At, = tj - tj-i: 
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10 



QgT = A;; 



a 



g- 



m 



(Q23)^/M^Bo) 



(42) 



Because <Q23^ and <Bo> can be determined by averaging over many baseline 
measurements, the noise contribution of the correction term in Eqn. 42 can be made 
much smaller than in the single baseline correction procedure of § 8.2. Since the 
weighting function of ag. is small compared to unity, the measurement noise in Eqn. 42 
will then be determined primarily by and a^., and will be approximately the same as 
for a trapezoidal filter of the same peaking time and about a factor of sqrt(2) better than 
by using Eqn. 38. Furthermore, because a baseline cMOS filter capture can now be made 
at any time the interval between successive steps is adequately long, rather than requiring 
one in conjunction with each charge cMOS filter capture, the method's deadtime is cut in 
half: to that of a single cMOS filter captiu-e. 



15 



20 



25 



8.4. Measuring the non-ideal terms O22 and Bn 

An effective way to measure Q23 at time ti is by making a baseline cMOS 
capture, per Eqn. 36, with the difference now that the times ti need no longer be 
correlated with charge cMOS captures. Using an averaged value <Bo> for Bo, plus 
computed values for L and M, we can then solve [LQ23 + MBojb for Q23: 



Q23 =1. 



Mr 



Jb 



(43) 



and use it to update <Q23>i using Eqn. 40. 



We normally expect Bo to be stable and nominally rate independent, so 
that <Bo> can be measured with no input coimts to the detector and used thereafter to 
obtain Q23 values. However, to compensate for any temperature and/or rate related 
changes in Bq, it is prudent to also measure Bo from time to time and update our estimate 
of <Bo>. One approach is to use pairs of baseline measurements, separated in time by 
At, where, if we denote the error term E as: 



E = 



LQ23+MBQ 



Jb 



A. 

A 



Jb. 



(44) 



then, measuring Ej and Ei.i separated by At = ti - ti-i, we have, from Eqn. 39: 
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Bo,= (E,-b^E,_„)(M(.-.f.))- ^^^^ 

Bo,i measurements can be averaged, as noted above, to obtain <Bo>, a low variance 
estimation of Bq, for use in Eqn. 35a. From Eqn. 45, it is clear that the longer the time 
interval At, the more accurate the measurement of Bq will be, since both the numerator 
5 and denominator are small differences of large numbers. In any case, enough values Bo,i 
can be averaged so that errors in <Bo> do not contribute substantially to the error in QgT- 

8.5, Dealing with model parameter drifts 

For this compensation scheme to work well, the P4 parameter p needs to 
be precisely known and stable over time. Else both the extrapolations of Q13 and 
10 increments in Q23 will incorporate systematic errors which will propagate into the Qgj 
calculations. The proposed Q23 baseline measurements offer a method to monitor p and 
either refine its value or track it in time. 

In Eqn. 40 we update our average <Q23>i = b2^^<Q23>i-l at time tj with 
the measured value Q23(i). However, we can also record the differences between these 
15 two values and compute their mean <AQ23> over time. If <AQ23> is not zero, then a 
systematic error exists and we can adjust P using: 

5p = <AQ23>/Sp, (46) 

where the slope Sp = d<AQ23>/dp is previously obtained from calibration mns, made at 
system setup time, wherein p is systematically adjusted and the resultant <AQ23> values 
20 recorded, so that Sp can be calculated. Similar methods can be used to monitor and 
correct other P4 parameters such as xi and T2, if appropriate. 
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10. Conclusion 

The foregoing description of specific embodiments has been presented for 
purposes of illustration and description. It is not intended to be exhaustive or to limit the 
invention to the precise forms described, and obviously, many modifications and 
5 variations are possible in light of the above teaching. These embodiments were chosen 
and described in order to best explain the principles of the invention and its practical 
application to thereby enable others in the art to best utilize the invention in various 
' embodiments and with such modifications as best suit the invention to the particular uses 
contemplated. 

10 As a first example, while all the described embodiments used trapezoidal 

or triangular filters, there is no such restriction in the technique. As shown, any filter 
describable by a set of weighting constants applied to successive signal samples, or its 
analog equivalent, can be used. The method will therefore work with essentially any 
realizable filter. Second, while the specific embodiments employed two, three, or four 

15 filters, there is no specific limit on the number of filters which can be used. Even a single 
filter can be used if its output is correctly sampled at several times. If the number of 
filters N exceeds the number of base amplitudes M, then matrix inversion methods can be 
used to find least squares best estimates of the M values, rather than exact values. Third, 
while all the filters used in the specific embodiments were symmetric, this is not a 

20 requirement of the method, as may be seen from the derivations. Fourth, while specific 
configurations of electronic circuits were shown in the implementations, other 
configurations can readily be devised which accomplish the same functions. The number 
of ADC bits or value lengths captured to the registers could be varied, for example. Fifth, 
even the use of electronic circuits is not required to implement the method. As per the 

25 initial derivations, for example, if successive ADC values are instead read into computer 
memory, all the subsequent processing steps can be carried out in software. Sixth, the 
method can also be applied to signals which are formally equivalent to outputs fi^om 
nominally single pole filters, whether the signal source is actually a single pole filter or 
not. Superconducting bolometers produce such output signals, for example. 

30 Therefore, the above description should not be taken as limiting the scope 

of the invention, as defined by the appended claims. 
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WHAT IS CLAIMFH T5;- 

1. A method for determining the integrated area Qgx of a pulse-like input 
signal by measuring a step-like output signal provided by a nominally smgle-pole (N-IP) 
device in response to said pulse-like input signal, the method comprising: 

applying a filter set having one or more filters to said N-IP device output; 

detecting the presence of a step-Uke feature in said output signal; 

in response to detecting said feature, capturing a set of correlated multiple 
output sample values (the area cMOS) fi-om one or more filters in said filter set; and 

forming a weighted sum of the sample values in said area cMOS to 
determine said integrated area Qgj (the determined area) of said input signal; 

wherein the weights in said sum (the area weights) are selected to 
compensate said detemiined area for errors arising either fi-om the time structure of said 
pulse-like input, or fi-om deviations in the N- IP device's response fi-om an ideal 
single-pole response, or fi-om both. 

2. The method of claim 1 wherein said pulse-like input signal has a finite 
time extent and a possibly time- varying amplitude as well. 

3. The method of claim 2 wherein:. 

both DC and long decay constant terms in said N-IP device's output signal 
are ignored in computing the weights used in determining the integrated area QgT ; 

a longer time constant trapezoidal filter is used primarily to estimate the 
amplitude of said step-like feature; 

an intermediate time constant trapezoidal filter is used to estimate the 
device's first-pole response to prior pulse-like inputs; and 

a shorter time constant triangular filter is used to measure the risetime 
characteristics of said step-like feature. 

4. The method of claim 2 wherein:. 

both DC and long decay constant terms in said N-IP device's output are 
ignored in computing the weights used in determining the integrated area Qgj ; 

two longer time constant running average filters are used primarily to 
estimate both the amplitude of said step-like feature and the device's first-pole response 
to prior pulse-like inputs; and 

a shorter time constant running average filter is used to measure the 
risetime characteristics of said step-like feature. 



32 



wo 01/17103 



PCT/USOO/23333 



1 5, The method of claim 1 wherein said N-IP device is not ideal, but has 

2 additional poles and/or zeroes in its transfer function and possibly a DC offset as well. 

1 6. The method of claim 1 wherein the number of individual sample values 

2 in said area cMOS exceeds the number of fihers in said filter set. 

1 7. The method of claim 1 wherein the number of individual sample values 

2 in said area cMOS is equal to the number of filters in said filter set. 

1 8. The method of claim 1 wherein the output of said N-IP device is 

2 digitized by an analog-to-digital converter, portions of the digitized output of which are 

3 transferred to a memory, and at least one of said appljdng a filter set, detecting the 

4 presence, capturing, and forming a weighted sum is carried out by a computer operating 

5 under software control. 

1 9. The method of claim 1 wherein said area weights are determined 

2 analj^ically by: 

3 developing a mathematical model of said N-IP device's output signal in 

4 response to said pulse-like input signal in terms of: 

5 i) a first (Pi) parameter which is the desired integrated area Qgx, 

6 ii) a set of one or more second (P2) parameters which describe said N-IP 

7 device's residual amplitude response to any previous pulse-like inputs and/or its response 

8 to the finite arrival time of the present pulse-like input, the number of said Pi and P2 

9 parameters being less than or equal to the number of sample values captured in said area 

10 CMOS, and 

1 1 iii) a set of one or more third (P3) parameters which describe the transfer 

12 functions of the said N-IP device and the filters in said filter set; 

13 describing by its weighting fimction each filter fi-om which a value is 

14 captured at least once in said area cMOS; 

15 applying said weighting function to said parameterized output model to 

16 produce a linear equation between each of the values captured in said area cMOS fi-om 

17 said filter and said Pi and P2 parameters; and 

18 solving the resultant set of linear equations to obtain coefficients for 

19 expressing said Pi and P2 parameters in terms of the sample values of said area cMOS. 

1 10. The method of claim 9 wherein, if a P3 parameter value is not well 

2 known, determining its value by: 
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3 selecting a set of values for the P3 parameter which bracket its expected 

4 value; 

5 for each member of said set of bracketing values, carrying out the method 

6 of claim 1 .6 to produce a set of coefficients expressing said P] and P2 parameters in 

7 terms of said captured sample values ; 

8 using said set of coefficients to collect data to measure the value of said Pi 

9 parameter and determine the uncertainty in said measurement; 

10 using said uncertainties determined for the members of said set of 

1 1 bracketing values to plot or otherwise parameterize said uncertainty in terms of said 

12 poorly known P3 parameter; and 

1 3 using said plot or parameterization to select the value of said poorly known 

14 P3 parameter which minimizes the measurement uncertainty in said Pi parameter. 

1 11. The method of claim 1 wherein, if capturing said area cMOS would 

2 violate causality for any fiher in said filter set, relative to detecting the presence of said 

3 step-like feature, increasing said filter's propagation delay sufficiently so that causality is 

4 not violated by said capturing. 

1 12. The method of claim 1 wherein timers are used to capture values fi-om 

2 said filter outputs as they appear so that capturing an area cMOS can be carried out 

3 asynchronously. 

1 13. The method of claim 1 wherein additional propagation delays are 

2 added to said filters so that capturing an area cMOS can be carried out synchronously. 

1 14. The method of claim 1 wherein said N-IP device output is digitized 

2 by an analog-to-digital converter; 

3 the filters in said fiher set are digital fihers; 

4 the outputs of said digital filters are captured into registers; and 

5 the determination of said integrated area Qgx is carried out by a digital 

6 computing device that reads the values captured in said registers, multiplies each such 

7 value by a weighting coefficient, and adds them together. 

1 15. The method of claim 14 wherein said digital filters are either 

2 triangular or trapezoidal filters. 

1 16. The method of claim 14 wherein:. 

2 a cMOS value aj;, captured fi-om a trapezoidal filter, is used primarily to 

3 estimate the integrated area Qgx; 
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4 a cMOS value 032, captured from a triangular filter, is used to estimate the 

5 N-IP device's first pole residual responses to previous pulse-like inputs; and 

6 a cMOS value ag, captured from a triangular filter, is used to probe the 

7 risetime region of said step-like feature; and 

8 said integrated area QgT is determined according to 

^ QgT + ^127^ ^\l^32 

1 17. The method of claim 14 wherein: 

2 a cMOS value as, captured from a trapezoidal filter, is used primarily to 

3 estimate the integrated area QgT; 

4 a cMOS value a 12, captured from a triangular filter, is used to estimate the 

5 N-IP device's first pole residual response to previous pulse-like inputs; 

6 a cMOS value a34, captured from a triangular filter, is used to estimate the 

7 N-IP device's second pole residual response to previous pulse-like inputs; and 

8 a cMOS value ag, captured from a triangular filter, is used to probe the 

9 risetime region of said step-like feature; and 

10 said integrated area QgT is determined according to 

I _j ^g I I 

QgT = JTiCT5:+ Jj2— +JT3^43"*' 

1 18. The method of claim 14 wherein one or more of said digital filters is a 

2 running average filter. 

1 19. The method of claim 14 wherein digital timers are used to measure the 

2 time intervals between said detecting a step-like feature and said capture of the individual 

3 sample values in the area cMOS. 

1 20. The method of claim 1 wherein the need for timing precision in 

2 capturing said area cMOS relative to the position of the risetime region of said step-like 

3 output signal is reduced by: 

4 designating a risetime probe filter whose time constant is significantly 

5 longer than the longest expected risetime; 

6 adjusting the area cMOS capture so that the risetime regions of the N-IP 

7 device's step-like outputs will fall entirely within it; and 

8 for all other filter values that are captured, adjusting the area cM OS 

9 capture so that the weighting constants their respective filters are negligible within said 
1 0 risetime region. 
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1 21 . The method of claim 20 wherein if one of said other filter values is 

2 captured from a trapezoidal filter, adjusting the area cMOS capture so that the risetime 

3 region lies within its gap and if one of said other filter values is captured from a triangular 

4 filter, adjusting the area cMOS capture so that the risetime region lies entirely outside its 

5 range. 

1 22. The method of claim 1 wherein the time of detecting said step-like 

2 feature is captured by: 

3 connecting a digital register to the output of a digital clock; and 

4 using said detecting the presence to cause said register to load the time 

5 value of said clock at that instant. 

1 23. The method of claim 1 wherein. 

2 the filters in said filter set are analog shaping filters; 

3 their output values, delayed if necessary, are captured by sample and hold 

4 circuits which are triggered by the outputs of timing circuits which are started by said 

5 detection event; 

6 said sample and hold outputs are weighted by constants and added using 

7 an analog weight-and-add circuit or multiply-and-add circuit; and 

8 if additional pileup criteria are satisfied, the output of this circuit is shaped 

9 by an output gating circuit so that its amplitude correctly represents the desired integrated 
10 area Qgx- 

1 24. The method of claim 1 wherein, in addition to detecting the presence 

2 of a step-like feature, an additional test is made to see if said pulse-like input signal is 

3 piled up and, if it is, preventing both the following capture and determination from being 

4 carried out, 

1 25. The method of claim 1, further comprising capturing the peak 

2 amplitudes of one or more filters and treating these captured values as members of the 

3 area cMOS in determining said integrated area Qgj. 

1 26. The method of claim 1, further comprising, when said N-IP device is 

2 not strictly a single-pole device, but has at least a second pole, plus possible additional 

3 poles and/or zeros and/or a DC offset: 

4 forming a mathematical model which describes the step-like output signals 

5 from said N- 1 P device in terms of 

6 i) a first (Pi) parameter that is the desired integrated area Qgx, 
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7 ii) a set of one or more second (P2) parameters that describe said N-IP 

8 device's single pole response of to any previous pulse-like inputs and/or its response to 

9 the finite arrival time of the present pulse, wherein the number of said Pi and P2 

10 parameters is less than or equal to the number of sample values captured in said area 

11 CMOS; 

^2 iii) a set of one or more third (P3) parameters that describe the non-ideal 

13 portion of said N-IP device's output, including Q23, its second pole response to any 

14 previous pulse-like inputs, and possibly Bq, its DC offset, as well (the error terms); 

15 iv) a set of one or more fixed fourth (P4) parameters that describe the 

16 transfer functions of said N-IP device and of the filters in said filter set; 

17 designating one or more filters in said filter set as the baseline filter set; 

1 8 capturing, from time to time, one or more output values fi-om said baseline 

1 9 filter set (a baseline sample); and 

20 using one or more baseline samples to update at least one type P3 

21 parameter. 

1 27. The method of claim 26, further comprising: 

2 using said type P3 parameters to correct said integrated area Qgj for the 

3 effects of said error terms. 

1 28. The method of claim 27 wherein said correction is accomplished by 

2 adding a weighted sum of said type P3 parameters to said integrated area Qgj. 

1 29 . The method of claim 28, when said set of P3 parameters includes only 

2 Q23 and Bq's correcting said integrated area Qgj comprises: 

3 using ^m^[LQ23 + MB(J ^^^^ weighted sum of said type P3 parameters, 

4 _ wherein Q23 is evaluated at a time fixed with respect to said area cMOS 

5 capture and , L, and M depend only upon type P4 parameters. 

1 30. The method of claim 26 wherein at least one type P3 parameter update 

2 is accomplished by: 

3 using said baseline sample to produce a new estimate of said type P3 

4 parameter; and 

5 forming an appropriate average between said new estimate and said 

6 parameter's prior value. 

1 31. The method of claim 30 wherein, in forming said appropriate average: 

2 the difference between said new estimate and said prior value is recorded 

3 and averaged as well; and 
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4 this difference average is used to adjust one of the type P4 parameters. 

1 32. The method of claim 26, further comprising: 

2 recording the time of capture of each baseUne sample (time ti); and 

3 using ti in said parameter updating to compensate for said P3 parameters' 

4 variation in time, so generating a time compensated parameter model (TC-model). 

1 33. The method of claim 32 wherein the update of an averaged estimate 

2 <Q> of a time varying type P3 TC-model parameter Q is accomplished by: 

3 using said baseline sample to produce a new estimate Qi of said parameter 

4 Q at said time ti ; and 

5 averaging Qi with said parameter's averaged estimate <Q>i.i at the prior 

6 measurement time ti.i to produce a value <Q>i of <Q> at time t\ according to: 

8 where N is a constant and b corrects <Q> for its time variation between times ti and ti.i. 

1 34. The method of claim 33 wherein the type parameter <Q> is associated 

2 with a pole having a decay time Xq and b is given by: 

3 l> = exp[-(t,-t,_,]/x^) 

1 35. The method of claim 32 comprising, whenever the integrated area Qgx 

2 of a pulse-like input is determined, using QgT to update said TC-model parameters to 

3 reflect changes in said error terms in response to said pulse-like input. 

1 36. The method of claim 35 wherein: 

2 a P3 parameter Q is associated with a pole having decay time Xq; 

3 said time of determining said integrated area QgT is denoted ti; 

4 the most recent prior time of updating Q, to the value Qi-i, in response to 

5 capturing a baseline sample is denoted ti.i; and 

6 said updating of Q to value Qi to reflect changes in said error terms in 

7 response to Qgj is accomplished according to: 

9 where the type P4 parameter p is the fraction of QgT which is transferred into the Xq decay 
10 mode. 
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1 37. In the method of claim 32, further comprising: 

2 recording, in said area cMOS capture, the time of said capture; 

3 adjusting said type P3 TC-model parameters for the elapsed time Ati 

4 between said time of area cMOS capture ti and the last time ti_i the TC-model was 

5 updated; and 

6 using said adjusted type P3 parameters to correct said integrated area Qgx 

7 for the effects of said error terms. 

1 38. The method of claim 37 wherein said integrated area correction step is 

2 carried out by: 

3 applying a set of weights to said type P3 TC-model parameters; and 

4 adding a sum of said weighted P3 parameters to said integrated area Qgj- 

1 39. The method of claim 38 wherein. 

2 said set of P3 parameters includes only Q23 and Bq; 

3 said area cMOS comprises three captured values, designated a£+, aj:-, 

4 and Og, with Og probing the step-like feature's risetime region; and 

5 said weighted sum of said area cMOS values used to determine said 

6 integrated area Qgj, minus said correcting weighted sum of type P3 parameters, is given 

7 by: 

QgT= -^m Lbf^(Q23)^^ ^(^o) 

8 ^ ^ ^ \ 

9 where 

10 K, , , L and M depend only upon type P4 parameters, 

1 1 <Bo> is the TC-model estimate of Bq, 

12 <Q23>i-l is the TC-model estimate of Q23 at its time of last update tf.i, and 

13 b2 is the amount Q23 exponentially decays per unit time. 

1 40. The method of claim 26 wherein one or more of the filters whose 

2 output values are captured in said area cMOS are also used in the baseline filter set. 

1 41 . The method of claim 26 wherein said baseline value captures are 

2 allowed only at instants when said baseline filter set is not engaged in processing a step- 

3 like feature in the N-IP device's output. 

1 42. The method of claim 26 wherein said baseline capture is a cMOS 

2 capture. 
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1 43. The method of claim 26 wherein: 

2 said N-IP device output is digitized by an analog to digital converter; 

3 the filters in said baseline filter set are digital filters; 

4 the baseline sample is captured by capturing the digital outputs of said 

5 filters into registers; and 

6 the model parameter update step is carried out by a digital computing 

7 device which reads the values captured in said registers. 

1 44. The method of claim 43 wherein said digital computing device carries 

2 out said correction of said integrated area Qgx by adding a weighted sum of said 

3 TC-model parameters to Qgj. 

1 45. The method of claim 26 wherein: 

2 said one or more filters fi-om which said area cMOS is captured is 

3 designated the area filter set; 

4 the said area and baseline filter sets are identical; 

5 identical correlated multiple output samples are captured fi-om said 

6 baseline filter set (a baseline cMOS) and said area filter set (i.e., the area cMOS); and, 

7 additionally, 

8 said baseline cMOS captures are allowed only at instants when said 

9 baseline filter set is not engaged in processing a step-like feature in the N-IP device's 

10 output; 

1 1 said area weights are also applied to the baseline cMOS values to compute 

12 null area estimates of the response of said area filter set to said second pole plus possible 

13 additional poles and/or zeros and/or DC offset; and 

14 one or more of said null area estimates are used in said updating of at least 

15 one type P3 parameter. 

1 46. The method of claim 45 wherein: 

2 said set of P3 parameters includes only Q23 and Bo; 

3 said area cMOS comprises three captured values, designated as+» crx-, 

4 and ag, with ag probing the step-like feature's risetime region, and 

5 said weighted sum of said area cMOS values used to determine said 

6 integrated area Qgj is given by: 

QgT = <^S+ - J^S- + ^^g- 

7 I ia; 
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values of Q23 are obtained from said null area measurement according to: 



Q23 



Mo 



Jb 



9 

10 where , K, , ^ , L and M depend only upon type P4 parameters, and 

1 1 the subscripts a and b stand for area cMOS captures and baseline cMOS captures, 

12 respectively. 

1 47. The method of claim 1, and further comprising, when said N-IP 

2 device is not strictly a single-pole device, but has at least a second pole, plus possible 

3 additional poles and/or zeros and/or a DC offset: 

4 designating one or more filters in said filter set as the baseline filter set; 

5 capturing, from time to time, one or more output values from said baseline 

6 filter set (a baseline sample); and 

7 using one or more of said baseline samples to correct said determined area 

8 QgT for errors arising Srom of that portion of the N-IP device's output (the error terms) 

9 . which is due to said second pole plus possible additional poles and/or zeros and/or DC 
10 offset. 

1 48. The method of claim 47 wherein: 

2 said one or more filters from which said area cMOS is captured is 

3 designated the area filter set; 

4 the said area and baseline filter sets are identical; 

5 identical correlated multiple output samples are captured from said 

6 baseline filter set (a baseline cMOS) and said area filter set (i.e., the area cMOS). 

1 49. The method of claim 48 wherein: - 

2 said baseline cMOS captures are allowed only at instants when said area 

3 filter set is not engaged in processing one of said step-like output signals from the N-1 P 

4 device; 

5 said area weights are also applied to the baseUne cMOS values to compute 

6 null area estimates of the response of said area filter set to said second pole plus possible 

7 additional poles and/or zeros and/or DC offset; and 

8 one or more of said null area estimates are used in said correcting of said 

9 determined area Qgj. 

1 50. The method of claim 49 wherein said correction of said determined 

2 area Qgj is carried out by: 
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optionally multipljong said null area estimate by a time decay factor b*^ 
where said coefficient b describes the expon^tial decay of said null area per unit time 
and 5t is the difference between the time of the area cMOS capture and the time of the 
baseline cMOS capture; and 

subtracting said null area estimate from said determined area Qgx- 

5 1 . The method of claim 50 wherein. 

said area cMOS comprises three captured values, designated CT2+, <7S-. 
and og, with CTg probing the step-like feature's risetime region, and 
said corrected area Qgx is given by: 



QgT=A-* 



As. 

A 



5t 



where the subscripts a and n label the area and null area estimates, respectively, and , 



A 



m 



K, and r are said area weights. 

52. The method of claim 49 wherein said correction of said determined 
area Qgx is carried out by subtracting an average of several null area estimates from said 
determined area. 



1 53. The method of claim 1 wherein compensation is made for a finite time 

2 extent of said pulse-like input signal by the additional steps of: 

3 in response to said detecting of a step-like feature, capturing a digital 

4 representation of the risetime portion of said detected step-like signal; and 

5 using a digital computing device to process the values in said 

6 representation in order to compute said compensation. 

1 54. The method of claim 1, 3, 4, 8, 12, 13, 14, 15, 18, 19, 20, 22, 23, 32, 

2 44, 45, 47 or 49 wherein said N-IP device is a charge sensitive preamplifier, said pulse- 

3 like input signal is a current from a radiation detector following a radiation absorption 

4 event, and said determined integrated area Qgj is a charge which represents the energy 

5 deposited in said detector by said absorption event. 

1 55. The method of claim 1, 8, 14, 15, 20 or 47 wherein said N-IP device 

2 is a superconducting bolometer, said pulse-like input signal is the flow of energy 

3 following absorption of a photon in the bolometer's active volume, and said determined 
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4 integrated area Qgj represents the energy deposited in said active volume by said photon 

5 absorption event. 

1 56. Apparatus for measuring a step-like output signal from a nominally 

2 single-pole device (N-IP device) in response to a pulse-like input signal to determine the 

3 integrated area Qgx of said pulse-like input signal, the apparatus comprising: 

4 a set of one or more filters applied to the N- IP device output (the area 

5 filter set); 

6 circuitry that detects the presence of a step-like signal in said output signal; 

7 circuitry that captures a correlated multiple output sample (the area cMOS) 

8 from said area filter set in response to said detection event; and 

9 circuitry that forms a weighted sum of the values in said area cMOS to 

10 detemiine said integrated area Qgx (the determined area) of said pulse-like input signal; 

1 1 wherein the weights in said sum (the area weights) are selected to 

12 compensate said area determination for errors arising either from the time structure of 

13 said pulse-like input signal, or from deviations in the N-IP device's response from an 

14 ideal single pole response, or from both. 

1 57. The apparatus of claim 56 wherein: 

2 said N-IP device output is digitized by an analog to digital converter; 

3 said filters in said area filter set are digital filters; 

4 said circuitry that captures include digital registers into which the digital. 

5 outputs of said filters are captured; and 

6 said circuitry that forms a weighted sum to determine said integrated area 

7 QgT includes a digital computing device that reads the values captured in said registers, 

8 multiplies each such value by an appropriate weighting constant, and adds them together 

9 to compute said integrated area Qgj . 

1 58. The apparatus of claim 56 wherein: 

2 said N-IP device is a charge sensitive preamplifier; 

3 said pulse-like input is from a radiation detector following a radiation 

4 absorption event; and 

5 said determined integrated area Qgx represents the energy deposited in said 

6 detector by said absorption event. 

1 59. The apparatus of claim 56, which further compensates for the finite 

2 time extent of said pulse-like input signal, the apparatus additionally comprising: 

3 circuitry that captures a digital representation of the risetime portion of 

4 said detected step-like output feature in response to said detection event; and 
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5 a digital computing device that processes the values in said representation 

6 in order to compute said compensation. 

1 60. The apparatus of claim 56 wherein, in order to produce a prescribed 

2 set of time relationships between the members of said area cMOS, said circuitry that 

3 captiu^es additionally includes, for each filter value to be captured, correlated sampling 

4 circuitry comprising one or more of: 

5 a timer that fixes the time of capture relative to said detection event; or 

6 a delay element inserted between the output of said filter and said circuitry 

7 that captures; or 

8 wherein the output of a single filter may be captured more than once by the 

9 use of multiple sets of said circuitry that captures. 

1 61 . The apparatus of claim 60 wherein: 

2 said N-IP device is a charge sensitive preamplifier; 

3 said pulse-like input is fi-om a radiation detector following a radiation 

4 absorption event; and 

5 said determined integrated area Qgj represents the energy deposited in said 

6 detector by said absorption event. 

1 62. The apparatus ofclaim 60 wherein: 

2 said N-IP device is a superconducting bolometer; 

3 said pulse-Uke input is the absorption of a photon in the bolometer's active 

4 volume; and 

5 said determined integrated area QgT represents the energy deposited in 

6 said active volxmie by said photon absorption event. 
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